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Assuming a tricritical point of the two-fiavor QCD in the space of temperature, baryon num- 
ber chemical potential and quark mass, we study the change of the associated soft mode along the 
critical line within the Ginzburg-Landau approach and the Nambu-Jona-Lasinio model. The or- 
dering density along the chiral critical line is the scalar density whereas a linear combination of the 
scalar, baryon number and energy densities becomes the proper ordering density along the critical 
line with finite quark masses. It is shown that the critical eigenmode shifts from the sigma-like 
. . . fiuctuation of the scalar density to a hydrodynamic mode at the tricritical point, where we have two 

\l , ordering densities, the scalar density and a linear combination of the baryon number and energy 

■ densities. We argue that appearance of the critical eigenmode with hydrodynamic character is a 

' logical consequence of divergent susceptibilities of the conserved densities. 
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I. INTRODUCTION 



O^l I At high temperature and/or baryon density, the system governed by QCD will show a transition from an ordinary 
. hadronic phase to a chirally symmetric, deconfined plasma phase [1,2]. The main objective of the heavy-ion programs 
at RHIC and at future LHC is to create this long-sought plasma state and to study collective properties of this 
, many-body assembly [3]. These two phases would have to be separated by a boundary with singularity if chiral 
' symmetry and/or confinement of QCD were exact symmetry. In reality, dynamical quarks with finite masses m 7^ 
, make both symmetries only approximate, and their order parameters, the quark condensate and the Polyakov loop, 
' have non-vanishing values everywhere in the phase diagram. Thus the plasma state may be smoothly connected with 
\ the ordinary hadronic state, even though they would possess qualitatively different properties from each other. 
(— I • Recently a strong possibility of a critical point in the real QCD phase diagram was suggested [4-6], based on model 
O 1' calculations [7-14] as well as lattice QCD results [15-17]. It is the endpoint of the first-order line, inferred from the 
JL : crossover behavior along the temperature (T) axis and the first order transition along the axis of the baryon-number 
chemical potential (/xs), and is a genuine singular point with the same criticality as the Z2 Ising model. Its location, 
(— I which is sensitive to the strange quark mass TOg, is expected to be within the reach of current experimental facilities. 
I!" ' Observable implications of this Z2 critical point (Z2CP)^ in heavy ion experiments have been discussed in the literature 
[4-6,18-23] such as large fluctuations of the low-momentum particle distributions, and the limitations on them due 
to the finite space-time geometry of collision events. This Z2CP will become a critical cornerstone in the QCD phase 
^ ' diagram once its location is confirmed in experiments. 

Based on the approximate chiral symmetry, the scalar density is usually taken as the order parameter of the 
Ginzburg-Landau (GL) effective potential to describe the critical behavior at the Z2CP. In this description all the 
singularities associated with the Z2CP seemingly originate from softening of the scalar density fluctuations as the 
effective potential becomes flat there. Especially, it might be concluded that the sigma meson becomes massless as 
an immediate consequence of this critical point. 

As a basic fact, however, we should strictly distinguish between the chiral critical point with m = and the Z2CP 
with m 7^ — even within the chiral effective models. When the chiral symmetry is exact, the T-jjLs plane is divided 
into two domains of the symmetric and broken phases with a boundary line. But the symmetry argument is unable 
to fix the order of the singularity of this line, especially the possible existence of the tri-critical point (TCP) on this 
line. Since the Z2CP at finite quark mass is the remnant of this TCP, the relation of the Z2CP to the chiral symmetry 
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is very obscure. In fact, Z2 is not the symmetry of the imderlying interactions, but of the thermodynamic potential 
at this particular isolated point in the T-fis phase diagram [24]. From this point of view the Z2CP is different from 
the cliiral critical point. 

The flat effective potential indicates the large fluctuation of the scalar density. At the chiral critical point, this is 
related to softening of the sigma meson mode, so as to form a symmetric excitation spectrum together with the pions. 
On the other hand, there is no such symmetry constraint at the Z2CP. Actually the scalar fluctuation linearly mixes 
with fluctuations of baryon number density and energy density, and therefore not only the sigma mode but also the 
hydrodynamic mode are to be taken into account there to study the associated soft mode. Consequently the fate of 
the sigma meson mode at the Z2CP is non trivial. 

Recent calculations of the dynamic mode in the scalar channel using the chiral models [13,25] indeed showed that 
the sigma meson is massive at the Z2CP. Furthermore, another scalar mode with space-like momentum dispersion is 
identified as the soft mode associated with the Z2CP in the Nambu- Jona-Lasinio (NJL) model [25]. In this paper we 
shall confirm the result of Ref. [25] on the more general ground using the time-dependent Ginzburg-Landau (TDGL) 
approach, and extend the study to discuss the changeover of the soft modes along the critical line in the T-/iB-m 
space within the TDGL approach as well as the NJL model. 

Our investigation is based on two fundamental observations about the Z2CP. The first point is that the proper 
ordering density at the Z2CP is a linear combination of the scalar, baryon number and internal energy densities 
[24-26] , as mentioned above. Because of this mixing all the susceptibilities of these densities diverge with the same 
critical exponent at the Z2CP. In contrast, in the chiral critical transition, the susceptibility of the scalar density 
diverges with exponent 7 of the 0(4) model in the two-flavor case, while the other susceptibilities of the baryon 
number and the energy have the smaller exponent a. 

The second is a consequence on the dynamics following from the conservation of the baryon number and the energy. 
The fluctuations of these conserved densities^ are intrinsically soft and constitute the hydrodynamic modes, whose 
excitation energies vanish as the wavevector q goes to zero. Susceptibilities of these conserved densities in turn 
have the spectral contributions solely from these hydrodynamic modes when expressed as a sum of mode spectra 
[27-29,26]. Hence the divergence of the susceptibility of a conserved density must be accompanied by critical slowing 
of a hydrodynamic mode. The spectral contribution from this hydrodynamic mode may well be involved in the scalar 
susceptibility through the mixing at the Z2CP. 

At an 0(4) critical point (0(4)CP) the importance of the hydrodynamic mode depends on which phase we start 
from. The hydrodynamic mode plays no critical role in the symmetric phase whereas the scalar condensate makes 
the mixing possible in the broken phase. The situation becomes more subtle at the TCP, where the 0(4) critical line 
shifts to the first-order line. Only the scalar susceptibility diverges due to the softening of the sigma meson at the 
TCP if it is approached from the symmetric phase. Otherwise, the hydrodynamic soft mode causes the divergence in 
the susceptibilities of the baryon number and energy as well as the scalar one. 

This paper is organized as follows. In the next section we briefly review generic properties of the phase diagram of 
QCD with two flavors near the TCP using the GL effective potential. It is stressed that at the TCP there are two 
relevant order parameters, the scalar condensate and a conserved density which is a linear combination of the baryon 
number and entropy densities. Then we include the dynamics using the TDGL model. Writing the susceptibilities as 
a spectral sum, we discuss the relative weight of the spectral contributions from the sigma and hydrodynamic modes. 
It is pointed out that the hydrodynamic contribution generates thc^ discontinuity of the baryon number and entropy 
susceptibilities at the 0(4)CP, and that this hydrodynamic mode gives the divergence at the TCP approached from 
the broken phase and also at the Z2CP. In §111 we perform the same analysis using the NJL model as an illustration. 
The GL effective potential with two ordering densities are numerically constructed there. The flat directions at the 
critical points are shown and discussed in relation to the divergences of the susceptibilities. The spectral origins of 
these divergences are studied with the relative weight of the mode spectra, and in detail based on the poles and 
residues of the scalar response function. Sections IV and V are devoted to discussions and summary. In Appendix A 
we prove the relation between the susceptibility and the response function, and in Appendix B we present the explicit 
formulas of the response functions in the NJL model. The results with the chiral quark model is briefly reported in 
Appendix C. 



^Momentum density is neglected here for simplicity. 
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II. GENERIC ANALYSIS 



A. Structure of the phase diagram and order parameters 

Let us briefly review the phase structure near the TCP [30,5,14]. It is known that the critical properties near the 
TCP are described, up to logarithmic corrections, with the Ginzburg-Landau effective potential 



ft 



= flo{T, Hb) + j d^x {a{T, i^b)(t^ + b{T, hb)ct^ + c{T, /xb)ct^ - ha) 

= QoiT, ijlb) + j (fxf{T, iJLB, h; a), (1) 

where / {Qq) denotes the (non-)singular part of the effective potential, and c > 0. The pseudo-scalar density is set 
to zero and neglected here in the mean field approximation. The critical exponents can be easily found from (1) at 
the mean field level. Along the line of the first-order transition within the symmetry plane (/i = 0), we have 

f = aa^ + ba^ + ca'' = ca\a^-alf, (2) 

where three minima with tr = 0, ±ai{T, iib) coexist (dashed line in Fig. 1). The baryon number and entropy densities 
are functions of cr^ due to symmetry, and discontinuous across the boundary between the symmetric phase (cr = 0) 
and the broken phase (cr = itri). At the TCP, where a = 6 = 0, these three phases coalesce and the first order line 
{h = ^2y^) smoothly joins with the 0(4) critical line (a = 0, 6 > 0). 

Once a small explicit breaking field —ha is exerted, the 0(4) critical line disappears and the TCP is lifted to the 
Z2CP. The line of Z2CP as a function of h is determined by the condition, /' = /" = /"' = (' = d/da), which is 
solved for a negative b with a = 36^/5c, b = —a'^/5 and a = sign(/i)(|/i|/16)^/^. Two lines of Z2CP with /i ^ form 
the edge of the wing-like surface of the first order transition in the a-h-h space, and these lines connect smoothly to 
the 0(4) critical line at the TCP. This wing structure is mapped into the physical phase space of T, /is and h ^ m 
(see Fig. 1). 

The slope of the first-order boundary can be related to the discontinuities of the densities across the boundary via 
the Clapeyron-Clausius relation [9,12], 

dT _ ApB dT Act dh _ ApB 



duB As ' dh As' d^B Act 



with baryon number density ps and entropy density s. The chiral broken phases with a = ±^ —a/2b coexist within 
the symmetry plane h = 0, and accordingly there is no gap in pB and s across this symmetry plane. Only the scalar 
density a bifurcates as the ordering density at the 0(4)CP approached from the symmetric phase. Its correlations 
with the "energy like" densities vanish (cts) = (apB) = because of the symmetry in a ^ —a in the symmetric 
phase. 

From the relation (3) we know that all of a, Pb and s generally have discontinuities across the wing because there 
is no reason for any of these slopes to vanish once h^ 0. Let us discuss the energy like and ordering densities around 
the Z2CP [24,31]. First we introduce the "temperature -like" field as a vector tangential to the coexistence boundary. 
Then the energy-like density is defined as the thermodynamic variable conjugate to this temperature-like field. This 
density has no discontinuity in the vicinity of the critical point. Since the boundary is two-dimensional, there are 
two independent temperature like fields and correspondingly two energy-like densities. Next the ordering density is 
defined as the density whose correlations with the energy like densities vanish at the critical point approached from 
the "symmetric" phase along the temperature-like direction. The conjugate field of this ordering density is no longer 
tangential to the coexistence boundary. There is a single ordering density at the Z2CP, which is in general a linear 
combination of a, ps and s. Since all the susceptibilities of these densities include the same singular fluctuation, they 
diverge with the same critical exponent at the Z2CP. 

The coexistence wing is squeezed to be one-dimensional at the TCP, where two lines of the Z2CP and the line of 
the first order transition with h = merge and smoothly connect to the single 0(4) critical line. Thus at the TCP we 
have only one energy- like density, which will be a linear combination of ps and s. Accordingly there are two ordering 
densities from dimensionality. The obvious one is the scalar density a related to the chiral symmetry and the other 
(fi is another linear combination of ps and s representing the Z2 symmetry of the potential at this particular point. 
It is sometimes useful to construct the effective potential with two ordering densities, a and which become soft at 
the TCP. 
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(a) 



(b) 



FIG. 1. (a) Schematic phase diagram axound the TCP in the a-b-h space. The three critical hnes are shown in bold lines 
which meet at the origin (TCP). The curve of three-phase coexistence is drawn in a bold dashed line which ends at TCP. The 
two-phase coexistence surface are hatched by thin dotted lines, (b) The counterpart in the physical T-fj,-^ space (NJL model). 



The same observation can be made by looking at the susceptibihties directly. There are three fields h, a and b in 
the eflfective potential (1). The singular parts of the corresponding susceptibilities form a 3-by-3 matrix = h, a, b), 

1 2a 40-3 \ 

X^J ^Xh\ 2a Aa^ 8a^ (4) 



with Xh the scalar susceptibility, 



4a^ Sa^ 16a^ / 



^'^ V 92/1 2a + 126ct2 + SOcct^ ' ' 



where a takes the value at the extremum of the potential. When the 0(4) CP is approached from the symmetric phase, 
Xaa = Xbb = 0. In the broken phase the situation is different. The singular part of (1) gives a finite contribution 
Xaa = 1/(3^) a-s the 0(4)CP is approached with a^ = —a/26 0, although the divergent susceptibility is still the 
scalar one alone. The singular contribution to Xaa eventually blows up at the TCP approached along the 0(4) critical 
line. In fact, the TCP may be understood as a usual critical point with the ordering density ip conjugate to a, sitting 
on the chiral phase boundary. When the TCP is approached from the broken phase with b = h = 0, the scalar 
ordering density vanishes slowly a^ = — a/3c, and Xh and Xaa diverge like l/\a\ and l/-\/[a[, respectively, while Xb 
is still non-singular. Note that the divergence of Xaa at the TCP indicates the infinities in the baryon number and 
energy^ susceptibilities, or equivalently the isothermal compressibility and specific heat, respectively. 

All these susceptibilities in (4) diverge at the Z2CP; where a 0. With the finite condensate a we can diagonalize 
this matrix of the susceptibilities, leaving only one singular susceptibility. The resulting eigenvalues are (0,0,(1 + 
4cr2 + I6a^)xh), with eigenvectors *(— 2cr, 1, 0), *(— 4(7^, 0, 1), *(1, 2a, 4(t^), respectively. For small a or h, we see that 
the ordering density is approximately a linear combination of the densities a and ip. 

The fluctuations of these two ordering densities a and ip become large near the TCP as explained above, and should 
be included as the soft degrees of freedom, especially when we discuss the dynamic aspects. We generalize the free 



^The fluctuation of the energy density is a linear combination of those of the baryon number and entropy densities. In this 
paper we sometimes use the energy susceptibility and the entropy susceptibility interchangeably. 
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energy so as to have two ordering densities, 

0, = j d^x ^aoCT^ + b^a^ + ca^ + 7crV + ^f"^ - ha- j(p^ + flQ. (6) 

Coupling between a and ip must respect the underlying chiral symmetry and the simplest coupling is a'^ip. A flat 
direction of this potential appears at a critical point in the a-(p plane. In the case of the 0(4)CP/TCP it is in the 
a direction reflecting the symmetry while the direction will become a linear combination of the two densities at the 
Z2CP. Eliminating the density (p by dVl/dp = 7(7^ + ^ — j =0, we recover the original form (1) of the free energy 
with a = ao + 7j and b = bo — ^7^, up to an analytic term. 



B. Dynamics 



We may introduce the dynamics to the system described by the free energy (6) phenomenologically [32]. We have 
seen in the previous subsection that there are two ordering densities conjugate to the fields h and a at the TCP, and 
that a linear combination of these two densities will become the relevant ordering density at the Z2CP for small h or 
a. We should include at least these two densities in order to describe the soft dynamics. Furthermore, it is known 
that (non-linear) mode-mode coupling between the fluctuations of the ordering densities and other (non-critical) 
hydrodynamic modes are important in general to describe the dynamics in the critical region [32], which is beyond 
the scope of this work [33,34]. We will see, however, that the coupled system of the two ordering densities in the mean 
field approximation yields already a non-trivial result [26] . 



1. Mixing between scalar and conserved densities 



Deviation of the densities from the absolute equilibrium gives rise to time evolution of the system. Here we assume 
simple phenomenological equations of motion for densities a and ip as 

L„{idt)a = -^, L^{idt)ip = (7) 

where L„{idt) and L^{idt) are the differential operators. Appropriate forms of and L^p are unknown in this 

description. But as a strong constraint we know that the operator L^p must be consistent with the conservation of the 
density p and describe the hydrodynamic motion. As a typical hydrodynamic evolution, we consider here the diffusion 
motion L^p{idt) = —dt/Xc^ with wavevector q. Note that the diffusion is time-irreversible. We assume propagating 
motion L„{idt) = df /V for the scalar density, identifying this mode as the sigma meson which degenerates with the 
propagating pions at the 0(4)CP/TCP. Other possible forms are considered below in this section. The transport 
coefficients F, A > are treated as constants here. 

For small deviations a ^ a + a and ip ^ ip + ip from the equilibrium values, we linearize these equations of motion 
with respect to a and (p to obtain 

L^{idt) + n^^ J [ip J 

where ilccr — d^i}/6a6a\eq, etc. The soft eigenmodes of the system are determined by the condition 

= 0, (9) 



+ r(x^' + 472(72 + Kq2) 27a 

27(7 ^rAq2 -iu + Aq2 



where Xh is the scalar susceptibility given in (5), and we introduced a term K(Vcr)2/2 in Q. The eigenmodes for small 
q are found as w = iwo, Wd with 



r 



Z^!^ = ^ = _^t1 (10) 
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where Xj is the susceptibihty of the density ip. 

The eigcnmodc lUq is oscihating while o^d has the diffusion-hke hydrodynamic character. The ujo vanishes at the 
0(4) CP, being the critical eigenmode. Although the hydrodynamic mode lu^ is an intrinsic soft mode of the system, 
it does not show the critical slowing there. When the TCP is approached from the symmetric phase, the situation is 
the same. On the other hand, at the TCP approached from the broken phase, both frequencies slow down, reflecting 
the divergence of the susceptibilities, which seems reflecting the existence of two independent ordering densities there. 

At the Z2CP the susceptibilities Xh and Xj diverge with the same exponent due to the linear mixing. The propagating 
LJo is a fast mode there due to the non-zero condensate a = (/i/16)^''^, whereas the hydrodynamic slow mode Wd/Aq^ = 
—ixj^ becomes the critical mode {xj^ 0) associated with the Z2CP. This result is similar to the level-crossing 
phenomenon where the mode coupling makes the lower energy mode lowered further. The explicit factor of Wd 
stemming from the hydrodynamic character results in the larger dynamic critical exponent ^; = 4 in the mean field 
level, which makes c^-vi apparently slower than the mode. In contrast, at the 0(4)CP, the linear mixing is banned 
by the underlying chiral symmetry. 



2. Susceptibility as a spectral sum 
Inverse of the differential operator (9) with the retarded boundary condition is the response function 

. . rW /^A^ + i -27^ ^ 

which characterizes the time-dependent response of these densities to the external fields, h and j, within the linear 
approximation. The susceptibility is obtain(;d in the limit of {uj = 0, q ^ 0). The response function is analytic in the 
upper complex-w plane, which fact allows us to express generally the susceptibility as a sum of the mode spectra [27]: 

X(0,q) = i /— Imx(a;,q), (12) 

where a ultra-violet regularization is understood if necessary. This expression shows that the divergence at a critical 
point should come from an infrared enhancement of the spectral function because the spectral function itself is usually 
integrable. 

Using this expression we can examine the relative weight of each mode contribution to the susceptibility. In our 
case, the oscillating and diffusion modes give spectral contributions as 

XH = lim - [ -ImxHi^, q) = Xh ( _i ' , , + , , ) , (13) 

and 

= q'So ^ / 77^™^^' ('^. q) = Xj (0 + 1) . (14) 

Here the first term in the bracket originates from the poles and the second from Wd- 

First, we note that only the diffusion like lu\i pole contributes to the susceptibility of ip. This is a robust result 
following from the conservation of the density ip. Existence of the current j such that di(p + \7 • j = dictates that the 
frequencies of the modes contributing to the (p susceptibility must vanish as q goes to zero. We can formally show 
that the spectrum of the (p response function behaves as limq^o IiiiXj(w, q)/w ex 6{ui) (see Appendix A). Conversely, 
we can state that softening of the hydrodynamic mode must accompany the divergence of Xj ■ 
Second, the ratio of the Ud spectral contribution to the total in the scalar susceptibility, 

47V2 

Xfc +47^2 ^ 

goes to unity at the TCP approached from the broken phase and at the Z2CP, which means that the leading divergence 
of the scalar susceptibility is also generated by the cod spectrum at these critical points. Even at the 0(4)CP approached 
from the broken phase the u>d spectrum gives a finite portion of the divergence < ii < 1 since Xh^ ~ ct^ ^ 0. This 
result can be understood by rewriting the scalar response function as 
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Xh{u>,ci) = xi°'(w,q)-— ^ ^ -, (16) 

1 - ^avXj (w, q)f^<7vXft (w, q) 

where xIj"-* (ti), q) = l/{La-{uj + ie) + fla-a) and xf'\^i^) = ^/{Lip{uj) + flipip). The denominator expresses the hnear 
mixing between the "bare" a and ip modes through the couphng fla-^p oc a. Even though the couphng becomes smaUer 
as the 0(4)CP is approached from the broken phase, softening of the mediating "bare" a propagator provides 1/cr^ 
factor, which results in the finite mixing of the Ud mode in the scalar channel. This is a simple example indicating 
the importance of the mode coupling near the critical point. 

In summary, the a and l/5 fluctuations mix and form two kinds of eigenmodes, ujo and lu^- We find that along the 
0(4) critical line approached from the broken phase, the critical eigenmode shifts from the sigma-meson like lUq to 
the diffusion like mode at the TCP. In contrast, when we approach the TCP from the symmetric phase, the scalar 
susceptibility Xh is given completely by the critical cOo spectrum without any mixing of the ip fluctuation. At the 
Z2CP the LOo mode becomes a fast mode whereas the whole divergence comes from the critical softening of the Wd 
spectrum with the hydrodynamic character. 

3. Cases with other types of motion 

In more microscopic NJL model calculation in the later section, the mode with the hydrodynamic character is 
provided as the mode of landau-damping type, contrary to the macroscopic analysis in the previous subsection, and 
we may change the time evolution operator accordingly as L^p{uj) ~ —iuj/\\/c^. More generally, the sigma mode may 
be oscillating or relaxing [Lc{w) = —uj^/T or —iui/T) while the possible hydrodynamic fluctuation can be the diffusion, 
landau-damping type, or sound motion (L^(w) = —iuj/\ci^,~iuj/\y/^, or — w^/Aq^). Furthermore, couplings with 
other hydrodynamic modes, if exist, are to be included to describe the correct dynamic behavior of the system [32 34]. 
We stress here that our result on the critical eigenmodes at the critical points are independent on these ambiguities. 
The only important fact is that the operator has the hydrodynamic character, and therefore intrinsically soft mode 
of the system. In any of these choices, we find the eigenmodes for small q^ as 

iaH = -(x^'+47V), L^{u:) = -xf- (17) 

Note that the typical mode frequencies of the diffusion, landau-damping type and sound-like dispersions vanish w ^ 
as q — > 0, and satisfy the spectral property following from the conservation law 

^"^ L,(a,)Vx-^ ""^^"^^ aslql-O. (18) 
Finally we note that these hydrodynamic modes drops out in the w limit of the response function: 

X(0+, 0) = hm^ - / -^^^^^lmx{uj' , 0) = lim - f ^ . u:'b{J)x = 0. (19) 



III. NAMBU-JONA LASINIO MODEL WITH A TRICRITICAL POINT 

As an definite illustration, we shall study the spectral contributions of the collective modes at critical points in 
the NJL model, and confirm that the result is consistent with the TDGL approach. We remark here that, unlike in 
the TDGL approach, there are no bare bosonic modes. The bosonic modes are dynamically generated through the 
interaction between the quarks and their softening causes the divergences at the critical points. 

A. Effective potential and susceptibilities 

We analyze the simplest version of the NJL model [35-37] £ = q{i^ — m)q + g[{qq)^ + {qij5T"'q)'^] in the mean field 
approximation {{qq) = CT=const, {qij5T°'q) = 7r=0). The thermodynamics is described by the effective potential [38], 

/(Pk 1 
-^[E - Tln(l - n+) - Tln(l - n_)] + -^{2gaf, (20) 
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FIG. 2. Effective potentials of the NJL model at three critical points, (a) 0(4) CP, (b) TCP and (c) Z2CP. In the upper 
panels shown are the potentials (20) measured from the minima as functions of a single ordering density a. The middle and 
lower panels are the contour plots of the potentials (24) with two ordering densities, (cr, p) and (a, s), respectively. 



where n± = (e^'^^^F^) + I)-\ E = v^U^Tl?, M = m- 2ga, and ly = 2NfNc = 2 • 2 • 3 = 12 with Nf and the 
numbers of flavor and color, respectively. Here /i is the quark chemical potential. The true thermodynamic state is 
determined by the extremum condition, dH/da = 0, and the corresponding grand potential is Q{T,iJ.,m). We define 
the model with the three-momentum cutoff A and with the coupling constant gA^ = 2.5 which allows the TCP. In 
the following, all the dimensionful quantities are expressed in the units of A. 
Expansion of the effective potential around a = with m = gives rise to 

-^[k - Tln(l - nl) - Tln(l - n° )] + -(- - J°){2gar + ^I^i^Q'^f + ■■■ 

= n{T,f,,0;0)/V + a{T,t,)i2gaf + b{T,t,){2gaf + ---, (21) 

where the superscript indicates the quantity evaluated in the massless limit. The first term is the non-singular part 
of the free energy in the GL description. The integrals 7° and J*^ are given in Appendix B. The TCP determined by 
a = b = appears at Tt/A = 0.20362 and /Xt/A = 0.49558. 

As explained in §11, it is useful to introduce the effective potential with another relevant ordering density besides 
the a in studying the behavior of the quark number susceptibility and specific heat near the TCP and the Z2CP. From 
the physical grand potential f2(T, /x, m), we can construct the Landau effective potential with two ordering densities 
p and a in the following way: first we introduce the free energy F{T, p, a) via 

F{T, p, a)/V = 0(T, fi, fh)/V + flp- ma, (22) 

where p, = il{T, p, a) and m = rh(T, p, a) are defined by inverting the functions 
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FIG. 3. Effective potentials of the NJL model as functions of a single ordering 
density p at three critical points, (a) 0(4)CP, (b) TCP and (c) Z2CP. The dashed 
line indicates the critical density in each case. 
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FIG. 4. X/JM (solid line) along the 

0(4) critical line approached from 
the broken phase. Xmm approached 
from the symmetric phase is shown 
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Then introducing new parameters fx and m, we define the Landau-type effective potential as 

f2(T, /z, m; p, a) /V = F{T, p, a) /V - up + ma = 0(T, /x, m)/V + {p. - ij)p - [m - m)a. (24) 

The extremum condition for the densities p and a yields p = p, and m = to, recovering the physical grand potential 
fl{T,p,m). Use of the entropy density s instead of the quark number p is straightforward. It is known that the 
effective potential constructed in this way must be convex and cannot be defined in the mixed phase. Fortunately in 
the NJL model we can bypass this difficulty by supplementing the unphysical grand potential ^{T, p, m) defined with 
the unstable solutions of the gap equation, which corresponds to the non convex part of the potential fi. 

One can easily show that the susceptibilities Xij = —yd^il/didj = T,p,m) are equal to the inverse of 
the curvature matrix of the GL effective potential at the extremum point. Prom this fact the divergence of the 
susceptibilities at a critical point is related to the appearance of a particular flat direction in the GL effective potential. 

In Fig. 2 we show the effective potential with two ordering densities at three critical points, (a) 0(4)CP with 
{T,p) = (0.3419,0.3), (b) TCP, and (c) Z2CP with {T,p,m) = (0.1498,0.5701,0.01) in the units of A. The critical 
instability at these points is usually discussed using the effective potential (20) with a single order parameter a based 
on chiral symmetry, which is shown in the upper panels. The flat curvature of this potential means the divergence of 
the scalar susceptibility. At the 0(4)GP it is clear from the potential with two ordering densities {a, p), or {a, s), that 
the (T axis is indeed the symmetry direction of th(^ system. The densities p and s depend on cr^ when calculated from 
(20) with TO = 0. This fact is seen here as a quadratic bending of the potential valley. Thus the fluctuation of these 
densities are weaker than that of a, and the susceptibilities of the quark number and the entropy have the smaller 
exponent a. 

At the Z2CP, on the other hand, the flat direction of the GL potential is not parallel to the a axis in the a-p 
and <T s planes. The proper flat direction is a linear combination of the three densities of a, p and s, and all the 
susceptibilities of them diverge with the same exponent at the Z2CP. 

It will be very instructive to introduce the GL function with single ordering density by eliminating a by 30. /da = 
in favor of p, as shown in Fig. 3. The curvature at the extremum coincides with the inverse of the quark number 
susceptibility. In case of the 0(4) CP the curvature does not vanish, implying the finite susceptibility Xt^n- It takes 
different values depending on from which side we approach the equilibrium value of p. Since the and a^ terms 
of the potential (21) disappear at the TCP and p changes with cr^ along the potential valley (sec Fig. 2 (b)), the p 
potential becomes flat on the side corresponding to the broken phase as seen in Fig. 3 (b) . This indicates the critical 
point for p. On the higher density side, in contrast, the curvature is non-vanishing. At the Z2CP, the potential Fig. 3 
(c) is essentially the same as the potential (20), and we may equally well choose p or s as the ordering density instead 
of a to describe this criticality. 

We show in Fig. 4 the quark number susceptibility as a function of p, along the 0(4) critical line, across which 
X/t/i is discontinuous. The value of Xnn on the 0(4) critical line approached from the broken phase grows up toward 
the TCP and eventually diverges there as is described with the GL potential (Xti/j. oc 1/6) [14]. The xtt also behaves 
in the same way. The only qualitative difference is that at the point p = the X/j-ti is continuous across the phase 
boundary because no linear coupling with a is allowed due to the symmetry under p ^ —p. 
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FIG. 5. Spectral functions of the scalar channel in the oj-q plane near the chiral critical point (Tc/A, //c/A)=(0.3419, 0.3). 
(a) T/A = 0.350 and (b) T/A = 0.339 with n = Axed. 



B. Response functions and mode spectra 



The spectral origin of the critical divergence can be investigated by studying the spectral function which is obtained 
as the imaginary part of the response function. We discuss here the structure of the collective eigenmode, which couples 
with the relevant susceptibilities and shows softening at the critical point. 

The response functions in the NJL model are calculated as [36-38] 

Xab{iqi, q) = n„6(%, q) + n„^(igi4, q)- ^ -2gYlmb{iqi, q), (a, b = fjL,m, (3). (25) 

Here the polarizations are defined with the imaginary-time quark propagator «S(^) = !/{/: + M) as 



^ab{iq4,cO = - j -0^Tj2tricBS{k)rsCk-q)r', 



(26) 



where (74 = 217tT {lei,), k — (k, k4 + i/x), F is an appropriate Dirac matrix, and the trace is taken over the flavor, 
color and Dirac indices. F = 1 for the scalar, 1^4 for the baryon number, and Hmf for with 

1 « 

Hmf = -i^liVi +M + (27) 

We calculate here the response function with Ti-MF the energy operator in the mean field approximation instead of 
the entropy, because the entropy has no microscopic expression as it is (lcfinc;(l only in equilibrium. The real-time 
response function is obtained from the imaginary-time propagator through the usual replacement iq4 ^ qo + ie in the 
final expression. The static response function in the long wavelength limit reduces to the corresponding susceptibility, 
limq^o Xafc(0, q) = Xab- Especially, xpp = T'^X.tt- 

The one-loop polarizations na6(a;, q) are the response functions of the free quark gas with mass M, and contain 
no contributions from the collective mode. In our simple NJL model, the collective mode is generated by the bubble 
sum encoded in the denominator of Eq. (24) and there are two kinds of collective motion [25]: the sigma meson mode 
and the particle-hole (p-h) mode. In Ref. [25] it is argued that the soft mode associated with the Z2CP is not the 
sigma meson but the p-h motion. 

The spectral function Pmm{^-,^ of the scalar response function [39] yields 

1 / 1 

Pmm{uJ, q) = 21mxmm(w, q) = 2Im— — - 1 

2.9 VI - 25(II„TO(a;,q} 



2Imn„„(w,q) 



[1 - 25Ren„„(a;, q)]2 + [2fflmn„„(a;, q)]2 ' 



(28) 
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FIG. 6. Spectral functions of the scalar channel in the w-g plane near the TCP. (a) T/A = 0.210 and (b) T/A = 0.2035 with 
fj, = Ht fixed. 

We notice that the spectrum 2ImnTOm('^) q) of the free quark gas is enhanced by the bubble-type correlation in the 
denominator. 

The scalar spectral functions arc shown at T/A=0.350 and 0.339 near the 0(4)CP with (Tc/A, ,Uc/A) = (0.3419, 0.3) 
in Fig. 5. One should keep in mind that the (5(w — 2M) spectrum of the sigma meson at q = in the broken phase 
is hard to be seen in this figure. The sigma meson spectrum is softening just above and below the 0(4)CP (see also 
§111 D). Besides the sigma spectrum we clearly find the p-h mode spectrum in the space-like momentum region, 
whose strength looks stronger in the broken phase. As we approach the TCP, as is shown in Fig. 6, this p-h mode 
spectrum grows in the small q region in the broken phase, while in the symmetric phase it does not show such a 
enhancement. 

We show in Fig. 7 the spectral functions of the scalar channel as well as the vector channel (quark number response) 

at the Z2CP with m,/A = 0.01. In the scalar channel clearly seen arc the two spectral peaks of the sigma meson 
and the p-h motion in medium, respectively [25]. This spectral structure is to be compared with that of the free 
quark case given in Appendix B. The most significant feature in Pmm is the critical enhancement in the Ct^ ~ region 
provided by the p-h mode, which gives rise to the divergence of the scalar susceptibility. Although the sigma meson 
shows the clear spectral peak in this model, the mode is massive due to the explicit symmetry breaking by the current 
quark mass. 

The spectral function of the quark number response, p^^(w,q), also contains these two spectral contributions, but 
the sigma spectrum strongly diminishes as q ^ 0. It is worthwhile to note the fact that limf^ioJ > 0, 0) = and 
Imnmp((x-', 0) oc luS{ll!) as q ^ 0, which reflects the conservation of the rjuark munber. Thus the response function 
X/*/i('^! 1) obtained in (25) in the random phase approximation (RPA), too, shares the same property, and the sigma 
meson cannot couple to the quark number susceptibility Xnn at q = 0. Therefore the divergence of Xij.h at the Z2CP 
must comes solely from the softening of the p-h spectrum. The same is true for Xtt- 

C. Spectral sum along the critical line 

In the previous subsection we have identified the soft mode associated with Z2CP as the p-h motion [25] generated 

in the scalar channel whereas at the 0(4) CP the sigma meson mode becomes soft. Let us examine how the changeover 
of the soft mode from the sigma meson to the p-h one occurs along the critical line. Once we noticed that the difference 
between the two limits, Xmm{0~^, 0) and Xmm{0, 0+), is caused by the hydrodynamic mode spectrum u}S{uj) as q ^ 0, 
it is easy to calculate the ratio of the spectral strength of the two types of spectral contributions. In these limits the 
explicit form of the RPA scalar response functions in the broken phase (with m = 0) yields, respectively. 
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FIG. 7. Spectral functions of the scalar (a) and quark number (b) response functions at the Z2CP with m/A = 0.01. 



Xmm(0+,O) = -^(^ 



2£f4M27(0+,0) 



1 



where function I{uj, q) is given in Appendix B. Then we can define the ratio R of the hydrodynamic spectrum to the 
total strength of the scalar susceptibility in the NJL model as 

„ ^ Xmm(0,O+)-Xmm(0+,O) ^ /(0+, 0) - J(0, + ) 

Xmm{0,O+) /(0+,0) • ^ ' 

On the other hand, these limits in the symmetric phase result in the same value 

which means no hydrodynamic contribution to the scalar susceptibility there {R = 0). The p-h mode must be 
decoupled from Xmm in the symmetric phase. Meanwhile we know that the corresponding ratios for the susceptibilities 
of the conserved quantities are always unity {R = 1), which can be explicitly seen with the expressions given in 
Appendix B. 

The ratio R (30) is shown in Fig. 8 as a fimction of /i along the critical line. We find that even in the 0(4) chiral 
transition at zero baryon number density (/U = 0) the hydrodynamic spectrum contributes to the divergence by a 
finite fraction. This contribution of the hydrodynamic spectrum increases toward the TCP, and eventually gives the 
leading divergence at the TCP, where 7(0,0"'") = 7'^ = but 7(0"*", 0) ^ 0. This behavior is completely in parallel 
with the TDGL approach. 

The fact that the p-h mode gives a finite fraction of the divergence at the 0(4) CP might be again unexpected from 
the viewpoint of the sigma meson as the associated soft mode there. Indeed, the sigma meson spectrum generates 
the total divergence when the critical point is approached from the symmetric phase. We should note here that 
the mixing of the scalar fluctuation in the broken phase is the origin of the discontinuity of the baryon number and 
energy susceptibilities across the boundary and that only the scalar p-h mode with the hydrodynamic character can 
couple with the fluctuations of these conserved quantities. Since the transitions between the scalar and other channels 
is proportional to M, the p-h spectral strength in the scalar is necessarily of order 1/M^ so as to bring a finite 
contribution to x^l^l■ xtt- 

We note that the scalar p-h motion of the NJL model is possible only in medium, but always possible in medium 
even in the symmetric phase, where the p-h contribution should be decoupled from the scalar susceptibility. One may 
ask the reason for this decoupling of the p-h spectrum. The absorption amplitude of the collective p-h mode with 
momentum q by a left-handed quark gi,(k) is proportional to a spinor product Ui{(k + q)ui,(k). In the symmetric 
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FIG. 8. Ratio (30) of the spectral contribution along the chiral critical line approached from the broken phase. 7i — > 1 toward 
the TCP. 



phase this coupling must contain a hehcity flip because chirality and helicity are identical for a massless quark, and 
vanishes in the q ^ limit. It is easy to show that UR(k + q)uL(k) oc |q|. For a massive quark in the broken phase, 
each of the left- and right-handed states has both helicity components, and thus uni^ + q)ML(k) 2M ^ in the 
q — > limit. 



D. Behavior of poles and residues 



Let us discuss a little more details of the spectral contributions to the susceptibility, studying the the poles and the 
residues of the scalar response function 

XmmM = ^ (i_25J(q) + 25(4M2-g2)7(^,q) ' ^ ^^^^ 
near the NJL critical points. It is useful to represent the spectral contributions as 

Z=±(T,ph 

with the poles corresponding to the sigma meson (iu) near ±2M and the p-h (ph) mode on the negative imaginary 
axis [25]. We would obtain the susceptibility as q — > after setting a; = 0. As for the p-h contribution, however, we 
take into account the kinematic condition \co/q\ < 1 for the spectrum via 

J_ ^'•'1 ^2Im ^P^^q) =^.l tan- ^. (34) 
27r i_|q| u> -UJ + Wph(q) ujph tt |a;ph| 

The scalar susceptibility in this approximation is expressed as a sum of the sigma and p-h pole contributions: 

X5?,t(0,0+) = 2^ + hm ^ . ^tan-1 (35) 

U>a q-*0 Wph TT |Wph| 

In Fig. 9 wc show numerical results of the polos and residues of the scalar response function (32) as functions of 
t = \T — Tc\/Tc with fixed = fic- The behavior of them can be understood as follows: 
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(T-TJ/T, (T-TJ/T, (T-TJ/T, 

FIG. 9. Scaling of the pole positions and the residues near the critical points, (a) 0(4)CP: The solid line denotes |a;ph|/|q| 
for i < and |wph|^/|q| for i > with q = 0.01. The residue \Rph\~^ for t < and \Rph\~^^^ for t > is shown in a dashed 
line. The sigma pole o;^ — 4M^ with q = for t < and (Rcci;^)'^ for f > in a dotted line, (b) TCP: The similar plot to 
the 0(4) case, but (|i^ph|/|q|)^ and for t < 0. (c) Z2CP: (|i^ph|/|q|)'^^^ x 10, \Rph\ and KeuJa are shown in solid, dashed and 
dotted lines, respectively. The IReJial is also drawn in a dash-dotted line. 



Across the 0(4)CP In the broken phase (1 — 2g J(0) = 0), the sigma pole with q = locates at lo„ = 2M on the real 
axis, whose residue is R^, ^ 1/(M/(2M, 0)). These quantities scale as \/Ra ^ M ^ Xmm ^ t^l"^. In the symmetric 
phase (M = 0), x^m 1 — 2gJ(0) ^ t. The complex sigma meson pole appears at u)„ ^ \J Xrnmli{^<y, 0) ^ t}^"^ 



with the residue ~ Xj ^Xmm.i{!^a^^) ~ i In both cases, the sigma mode gives the appropriate strength of the 

divergence Rjw^ ~ Xmm ~ t"^ 

The p-h mode arises from the w dependence of the function /(w, q). Since we are interested in the behavior in the 
small a; and q region with u = a;/|q| < 1, we may approximate the function I as 

/(w, q) = 7(0, 0+) + iulm7(0)', (36) 

where Im/(w)' = (a'/fiii)Im/(?i|q|, q)|q=o- In the broken phase, the condition, /(a),q) = gives the pole as tjph ~ 
— i|q|/(0, 0+), whose residue i?ph ~ — i|q|/M^ ~ Similarly in the symmetric phase the condition l/{2g) — 

J(0) — w^/|q| • iIm/(0)' = fixes the pole position as Wph ^ ~*(Xmml'll)^^^ ^ (^i^i)^^^ with the residue R^h ~ 
~*|Q|/(XmTOl*ll)^^^ ~ i~^^^|q|^^^- Then according to Eq. (34) the spectral contribution is estimated to be 

^tan-ii~i, (37) 
Wph |wph| t 

tan ^ - tan — fr- > (38) 



for the broken and symmetric cases, respectively, as |q| 0. We note that the p-h mode gives finite portion of the 
divergence at the 0(4)CP approached from the broken phase due to the enhancement of the residue by 1/M^, despite 
that the frequency Wph ^ |q| shows no critical slowing. The decoupling of the p-h mode in the symmetric phase is 
correctly described by the behavior of the pole. 

Across the TCP At the TCP, /(0+,0) ^ 1(0,0+) = 7° = 0. When we approach from the broken phase XmL = 
(2,9)24Af2/(0,0+) - t while ~ M ~ t^/*. Then the sigma mode slows down, but cannot generate the leading 
divergence because Ra/ijJa ~ t~^/'^ whereas Xmm ~ t~^- The fact 7(0, 0+) ~ t^/'^ changes the scaling of the p-h mode 
into Wph ~ t^^^|q| with the residue t~^/^|q|, which gives rise to the correct order of the divergence 

I^P^"tan-iM.ltan-^J_.i. (39) 



|wph| Wph t iV2 t' 

The p-h mode must correspond to the critical eigenmode. On the other hand, if the TCP is approached from the 
symmetric phase, the p-h mode is decoupled from the scalar susceptibility and the critical sigma mode generates the 
total divergence. It is very interesting that the soft mode associated with the TCP is different between the symmetric 
and broken phases. 

At the Z2 CP The sigma meson mode has a finite energy gap of order 2M in our model. The pole position of the 
p-h mode can be evaluated as Wph ^ — iXTOml^l ~ with its residue <^ — i|q|, which gives rise to 
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t^tan-iM^^_^f-2/3. (40) 

From this estimate, we see that the p-h mode properly accounts the divergence of the scalar susceptibility at the 
Z2CP, and therefore the softening of the p-h mode is the origin of the critical divergence at the Z2CP. 



IV. DISCUSSIONS 



In the microscopic calculation with the NJL model, we have seen that the collective p-h mode, besides the sigma 
meson mode, is generated in the scalar channel and brings the spectral contributions to the channels of the conserved 
quantities through the mixing when M ^ 0. This p-h spectrum makes the susceptibilities of and xtt discontin- 
uous across the 0(4) critical line, and eventually gives rise to the critical divergences at the TCP and the Z2CP. This 
role of the p-h mode is consistent with the behavior of the hydrodynamic mode in the TDGL analysis. We remark 
here that this p-h mode in the NJL model is the time-reversible landau-damping type. In the phenomenological 
TDGL approach, on the other hand, we assumed the time— irreversible diffusion motion for the conserved density, 
which seems more appropriate to the non equilibrium soft dynamics. It would be very interesting to study how a 
time-irreversible equation of motion emerges out of the time-reversible microscopic theory (see, e.g., Refs. [40-43]). 

The flat curvature of the effective potential is usually referred to as the vanishing screening mass, which naively 
hints the reduction of a kind of particle mass. As we approach the 0(4) CP, the sigma meson mass actually gets 
reduced to cause the critical divergence. However, approaching the Z2CP, we see that the flat curvature leads to the 
vanishing diffusion constant of the hydrodynamic mode [25,26,33]. In general, the potential curvature expresses the 
stiffness of the system with respect to the variation of the ordering density. The dynamic quantity related to this 
stiffness can be the particle mass, sound velocity, relaxation constant or diffusion constant, depending on the equation 
of motion of the critical eigenmode. 

In the Z2CP case, the linear mixing of the conserved densities in the proper ordering density dictates that the 
critical eigenmode should have hydrodynamic character. We have seen that the sigma-meson like mode is massive 
and is decoupled from the slow dynamics. It is explicitly argued in Ref. [33] that the remaining set of slow modes is 
equivalent to that of the liquid-gas critical point [25,26]. In the course of heavy ion events passing by the Z2CP, it is 
important to study the observable implications of critical slowing of hydrodynamic fluctuations in baryon number and 
entropy densities [5,18-23]. The growth of the diffusive fluctuations within the flnite space-time would get renewed 
interest [5,18]. 

The critical soft mode of the TCP is different between the symmetric and broken phases. The sigma mode becomes 

the critical eigenmode and the hydrodynamic one behaves normally if the critical point is approached from the exactly 
symmetric phase. Otherwise the critical eigenmode bears the hydrodynamic character and the sigma mode slows down 
only moderately, which implies the importance of the hydrodynamic fluctuations near the TCP: the endpoint of the 
flrst order line sitting on the chiral phase boundary. Theoretically and also experimentally it is worthwhile to elaborate 
and classify the dynamic critical behavior at the TCP as well as other critical points of QCD [26,33,34,44,45]. 

In the QCD thermodynamics the first order transition is believed to occur at finite temperature in the massless 
three-flavor case [2] and in the pure gluonic case [1], where the chiral symmetry and the center symmetry are 
exact, respectively. As varying masses of the quarks from these two limits, we will have critical points (line) at 
the edge of the first order region. For example, in Ref. [24], a Z2CP is studied in the T m. plane with three quark 
flavors of equal mass in lattice QCD. There a particular linear combination of the quark condensate and the energy 
density is identifled as a proper order parameter, which should be mapped to the magnetization in the Z2 Ising model. 
We expect that a hydrodynamic mode related to the energy fluctuation shows critical slowing at this Z2CP, whose 
spectrum may be detected in the lattice QCD. The dynamic critical behavior of this point is also of importance. It 
could be different from that of the Z2CP at flnite ijlb because of no linear mixing with the baryon number fluctuation 
due to the symmetry pB —ps- Furthermore, the lines and surfaces of the QCD critical points in the T-p,B-mud-^s 
space are speculated [46] . One can also extend the space to the isospin channel [47] . Generally at such a critical point 
the proper ordering density becomes a linear combination involving conserved densities. Since the critical eigenmode 
must have the hydrodynamic character in this case, the dynamic critical nature would be quite different from the case 
with (e.g.) the exact chiral symmetry. 

Our identification of the soft modes along the critical line is done within the mean field approximation. Fluctuations 
around the mean fields are known to become crucial for describing the singular behavior at the critical points correctly. 
To this end we should utilize mode-coupling theory or dynamic renormalization group method. The mean-field analysis 
provides a good starting point to identify an appropriate set of the slow modes. 



15 



V. SUMMARY 



The fundamental points about the Z2CP are following: (1) in the absence of the chiral symmetry, the ordering 
density becomes a linear combination of the scalar density, the baryon number density and the energy density, in 
general, and their susceptibilities have the same critical exponent. In describing the static property of this critical 
point one may equally well take any of these densities as the ordering density. (2) Then the critical eigenmode must 
be the hydrodynamic one which can cause the critical divergence of the susceptibility of the conserved density. On 
the other hand, in the chiral transition approached from the symmetric phase, the exact chiral symmetry prohibits 
the linear mixing of the hydrodynamic mode in the fluctuation of the ordering density. We have showed these points 
using the TDGL approach as well as the microscopic NJL model. 

We have studied the change of the critical eigenmode along the 0(4) critical line. When the critical point is 
approached from the symmetric phase, the soft mode is indeed the sigma meson mode. On the other hand, approaching 
from the broken phase we see the scalar condensate allows the linear mixing between the sigma and hydrodynamic 
modes, and eventually at the TCP the hydrodynamic mode turn out to be the critical eigenmode which generates 
the leading critical divergence. Thus the shift of the critical mode from the sigma meson to the hydrodynamic mode 
occurs at the TCP. And the soft mode at the TCP crucially depends on from which phase one is approaching the 
point. 

The criticality of the Z2CP is given by the softening of a hydrodynamic mode. The sigma mode remains as a fast 
mode due to the explicit breaking m and is decoupled from the soft dynamics. Based on this understanding, we should 
study fluctuations with the hydrodynamic character, such as baryon number and entropy fluctuations [5,18-23,48] in 
locating the Z2CP experimentally. 

In the QCD thermodynamics with three quark favors, several kinds of end points are speculated. One should keep 
in mind that at these points the hydrodynamic mode will become the critical eigenmode once a conserved density is 
linearly mixed in the ordering density and also that the approximate symmetry of the underlying interactions does 
not provide any reason for singularity in the phase diagram. 
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APPENDIX A: SUSCEPTIBILITY AND RESPONSE FUNCTION 



The susceptibility of a (bosonic) density 4>a in a system described with the grand potential Q 
with K = H — a J d^x4>a is defined as 



rintr(exp(-/3i^)) 



1 

Vdadb' 



(Al) 



This susceptibility is obtained as g'-limit of the response function [27,28] because 



d 



e^"tr(<^fa(0,O)e-'=*^) 





(A2) 



where we used a formula, ^e^^^") = dse^^^'^^ s)A(a) .^^j^jj g, matrix-valued function A{a), and the time de- 

pendence of the operators are defined by ^{—ir) = e'^^(f){0)e~^^ . The imaginary-time correlation and the response 
function are introduced as Xabi~ ~ {^a{—iT,x)^h{Q,0))c and Xab{t,x) = i6{t){[^a{t,x),(j)b{0,0)])c, respectively. 
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Once wc establish the relation between the susceptibility and the response function, it is easy to express the suscep- 
tibility as an integral over the spectral density: 



Xa6 = Xa6(0,q^0) 



lim 



du) 21mxab(a;,q) 

277 LO 



(A3) 



Specifically, when (j)a is conserved, [K, 4>a] = 0, we can freely change the position of the operator (j)a in the trace 
and the susceptibility is directly related to the equal-time correlation function S via [27,28] 



= (3 J d^xSab{0,^)=f3]imJ ^5ab(a;,q) 



:/31im / 

1-^0 J 



doj 21mxab(t^, q) 
2tt 1 - e-f^" 



(A4) 



In the last equality we used the fluctuation-dissipation theorem, which relates the fluctuation S to the dissipative 

part Imx of the response function. Noting the spectral condition, sign(oj)lmx{uJ . q) > 0, we conclude that these two 
expressions for the susceptibility coincide if and only if limq^o Imx(a;, q) — Tvd{uj)ujx- Physically this is a consequence 
of the existence of the current ja such that dt(l)a + V • ja = 0. 



APPENDIX B: EXPLICIT EXPRESSIONS 



1. response functions 



We present the explicit formulas and procedures to evaluate for the one-loop polarization functions. 
First working in the imaginary time formalism [38], we derive the expressions for the polarization functions with 
make use of the frequency sum formulas: 
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where = 21itT is the bosonic Matsubara frequency, fc = (k, - 



-L0„ 



- ifi) the quark momentum with w„ = (2n + l)7rT, 



fc| = k2 + kl El = x/M^n?, E2 = + (k 



q)2, and n±i^2 = n±{Ei^2)- The q- and k-dependences are 
implicit through these quantities. Then the analytic continuation to the real frequency is done by the replacement 
iuji ^ qo + is, which uniquely gives the retarded functions with the asymptotic behavior oc l/qo as l^ol ~* 00. 

Since the retarded function is analytic in the upper half plane, we can reconstruct it from the imaginary part using 
the dispersion integral. 



„ , , i r , , Im^„fc(^^q) 
n„b(w, q) = - / duj — r^. 



(B3) 



This relation is quite useful when we evaluate the polarization functions with finite q because the imaginary part is 
easier to calculate. The imaginary part comes from two physical processes in our model: the q-q creation/ annihilation 
and the mode-absorption/emission by a quark q or an anti-quark q as shown in Fig. 10. Kinematically the former 
occurs for the time-like momentum with = q^ — > AM"^. The latter is possible for q^ < 0, resulting in the 
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-2E -q q 2E Re 
X ^ ><p-h ^ ^ 



(a) (b) 

FIG. 10. (a) Pair creation/annihilation (upper) and absorption/emission (lower) processes contributing to the mode spectrum 
with the detailed balance, (b) Schematic analytic structure of the response function at finite momentum q. The cut (bold line) 
on the real axis with \oj\ > 2E is due to the upper processes of (a) while the cut between oj = —q and q is from the lower 
processes of (a) . The a meson pole and the p-h pole on the unphysical Riemann sheet are shown with x . 



Landau damping. Then the complex lo plane has two kinds of cuts in the case of finite q. Analytic continuation of 
the retarded function to the lower half plane (Im^; < 0) across one of these cuts is defined by 

n(z, q) = - r ^^ I™"(^'q) + 2Imn*'"'^P(0, q), (B4) 

TT W - Z 

where n*'"'''P(z, q) is the imaginary part in the time-like or the space-like region, depending on across which cut the 
function is continued. Substituting these expressions to Eq. (25), we can obtain the response functions with finite q 
and search the poles on the unphysical Riemann sheet (see Fig. 10). 
The one-loop polarization function in the scalar channel yields 

/ N f <Pk 1 — n+1 — n_i f (Pk , 9 , , ,ox^/ . x 
n„„(go, q) = ^ \ + V J j^{q^ - AM^)I{q, + ze) 

= J(|q|) + (g2-4M2)/(go + i£,q), (B5) 

where we have shifted the momentum as k ^ k -I- q/2 and £'1^2 = \/ M'^ -f (k ± q/2)2, and then introduce the cutoff 
at fc^ax = — q^/4 for this new k. In the second line wc defined the functions J and /, whose massless limits appear 
in the expansion series of the NJL effective potential in terms of a, Eq.(21). The integral l'^ is seemingly divergent 
logarithmically in the infrared region, but is actually finite because 1 — — n° ^ as k ^ 0. 

After performing the angular integration using the delta function imposed by the on-shell condition, the imaginary 
part of the scalar polarization for qo > yields 

X ^ / N V f dkk q^ - AAP ^, 



V{qo, q^) = \ ^ ~ '^+' " " " ^i^' ^ ^ ^^^^ ^ ''"^^"^ (B7) 



n+i - n+2 + JT— 1 - 'T— 2, q <0 

where = qf^ — q^ and E'^ ~ K'P + k^. Due to the on-shcll condition of the imaginary part, the quark energies in n± 
arc set to be £1,2 = f ± ^JW~q^ for q^ > AM'^ and Ei^2 = \/E'^ -q^/4± ^ for q'^ < 0. The last /^-integration 
can be done analytically. The imaginary part for q'o < is obtained by Imn(— qoiQ) = "111111(901 q)- 
Similarly, other polarization functions and their imaginary parts for qo > are found and given below: 

/ \ f d^k 1 - 77+1 - n_i f cPk , 2 , \ /■n.ox 

n^^(go,q) = y J +y J j^{ql-^E')I{qo), (B8) 

v f dkk qn — AE^ ^, „, 



imn. 
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/(Pk 



Imn„„(go,q) = - 



2tt) 
AMv f dkk 



167r J |q| 



(BIO) 
(Bll) 



(27r) 



+^ / [((90 - ^E^)iE' + M^) - EW + (k • q)2) I{qo) + n{ql - 4.E^)IM] 



Imn^/3(go,q) = 7^ / 

IbTT J 



V [ dkk 



(gg-4i?^)(g^-gV4 + /x^) 

_ ^2/4 



(B12) 
(B13) 



TT / \ w /" 1 - n+1 - n_i 
nm/3(go,q) =Mi/ y —-3 ^ +Miy 



El 



-^[{q'^-AE^)I{q,)-2t,IM], 



167r 


J ^ 


q) = 






-/ 


q) = 




16^./ 



d'^fc 1 — n+i — n_i 



!7r)3 

(2^)3 

r dkk 



El 



^E^ - g2/4 
In these expressions we introduced a function 

^w(go,g^) = ±(n+i +n_2 - n-i - n+2) 



(B14) 

(B15) 



(B16) 
(B17) 

(B18) 



with the sign '+' for AM'^ < < 4(M2 + fc^^^) and '-' for q^ < 0. 

It is known that the response functions have a non analytic property at the origin of the LD-q plane when the 
p-h mode spectrum exists. For demonstration we show in Fig. 11 the real and imaginary parts of Ilrnmi'^T'^ with 
q/A = 0.1. An abrupt change of the real part is seen in the region w/A < 0.1, which is clearly caused by the p-h 
spectrum in accord with the dispersion relation. In the q ^ limit, the imaginary part becomes proportional to 
uj5{lo), which leads to the discontinuity of the real part of Yi-mm as mentioned in Eq. (19) and shown in Fig. 11. In the 
case of the massless quarks, on the other hand, the scalar channel does not couple to the p-h motion in the q ^ 
limit due to the chiral symmetry. Hence the real part is non-singular at the origin as shown in Fig. 12. 

In Appendix A it is shown that the spectral function for a conserved density fluctuation must be proportional to 
u5{lo) as q ^ in general. We confirm that this requirement is fulfilled by these one-loop polarizations, H^^, 
11/3/3, and n„/3. 



2. susceptibilities 



The susceptibilities in free quark gas with mass M are found as 
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2lmn„„/A" 



(c) 




CO/A 



One-loop scalar polarization function Hmm (at the Z2CP). (a) The real part l/{2g) — Rellmm and the imaginary 
m are shown in solid and dashed lines, respectively, with fixed q/A = 0.1. (b) The real part as a function of {u!,q). 



FIG. 11. 
part Imn„ 
(c) The imaginary part as a function of 




yiO) 



A.mfj, 



X-mT 



y(0) 



(0) 

Xtt 



AO) 



(27r)3 

(Pk M . , , . 
J — (n+-n_), 



1 , ]VP fl-n+-n. 
-(l-n+-n_)-^( 



+ n+ + n_ 



f d-^ i 



(27r)3 E 
Sk M 



)3 E 

d^k . . . , 
(2^("++"-)' 

S'k 



E-^x , E + n , 
— — ^a. H = — n_ 



(27r)3 
d^k 



E- n 
T 



E + fi 
T 



f d'l 



E ~ fj, E + jj, 
— 7^ — n-L. = — n 



(B19) 



These expressions coincide with the static one-loop polarizations in the q ^ limit. Through this limiting procedure 
we find that the terms containing n'j_ are related with the p-h spectrum parts in functions T and T^i and that all 
susceptibilities, except for the scalar, must accompany because of their hydrodynamic nature. 

We also notice the fact of no mixing of the scalar fluctuation with others in the massless quark gas M = 0. In the 
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case of /i = 0, the vector fluctuation docs not mix with others. Both of these originate from the symmetries. If we 
write down the GL effective potential, it must be invariant under a —a and/or p — p, respectively, and therefore 
any linear mixing with other densities is impossible. 



APPENDIX C: CHIRAL QUARK MODEL 

The chiral quark model can be used to perform the same analysis as in the NJL model: 

1^X1 = - ^m'^(t>l - + ha + q[i^ - g{a + h5raT^a)]q, (CI) 

where 0o = = tto, and < 0. The meson mode is introduced here as an elementary field with the kinetic 

term. Integrating out the quark field, we obtain the effective potential within the mean field approximation for a and 

TT as 

n^^{T,n-<j)/V = -ha+-rh'cT^ + -a'-v j ^ - Tln(l - n+) - Tln(l - n_)] (C2) 



with E = V + and M = ga. This potential is almost the same as that of the NJL model (20), and this model 
is expected to have the same phase structure. The subtle point is that the divergent vacuum quark fluctuation in 
the integrand of Eq. (C2) requires a rcgularization and rcnormalization. Instead of the three momentum cutoff used 
in our NJL model calculation, here we adopt a simple prescription following (e.g.) Ref. [13]; we assume that the 
renormalization is already done in the vacuum and discard the vacuum polarization term. Then the parameters are 
chosen so as to reproduce the pion decay constant, the pion and sigma masses, and the constituent quark mass in the 
vacuum. We found the ZaCP at (Tc, Mc)=(117.7, 176.2) MeV. 

Within the same level of the approximation, the scalar response function is calculated as 

= -q^ + ni^ + kXa^-gm^^J ^^^^ 

where the polarization Hmm is defined in Eq. (26) with F = 1, but whose vacuum part is removed. The four-point 
interaction of the NJL model is replaced by the non-local one here. Other response functions 

Xab = ^ab + ^am „ , -0,1, o 577 ^mb {a,b = IJ,,0) (C4) 

-q^ +m^ + ^Xa^ - g^Tlmm 

have the same structure as the NJL result (25) because we assume the same scalar-type interaction between quarks. 

In the numerical calculation with this chiral quark model at the Z2CP, we confirmed the spectral enhancement in 
the space-like momentum region, similar to Fig. 7, and found a pole responsible for this enhancement on the negative 
imaginary axis in the complex-oj plane, just as in the NJL model. The ratio R defined in Eq. (30) also goes to unity 
as the Z2CP approached. Therefore our conclusion on the importance of the hydrodynamic mode at the Z2CP is 
unaltered here. 

We should note, however, that the semi-positivity condition on the spectral function is violated in our numerical 
result in the time-like momentum region. This is because we replaced the term 1 — n±i — nip2 in the expression of Tl-mm 
with —n±i — n^2 to remove the divergence (see Appendix B). This simple rcgularization breaks the detailed balance 
relation which is essential to assure thermal equilibrium. Hence the result of the spectrum in this chiral quark model 
should be interpreted with caution. The vacuum subtraction also results in the unexpected infra-red divergence of 
the quartic term, — A + 3g^I^ , in the expansion of (C2) around (7 = 0. Because of this difficulty, we could not find the 
TCP in this model with the rcgularization adopted here. In order to properly discuss the spectral structure in the 
chiral quark model we need the regularization scheme which satisfies the condition of thermal equilibrium [13,49]. 
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